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ABSTRACT  / 

The  popularity  of  micropropulsion  system  development  has  led  to  renewed  interest  in  the  determination  of 
propulsive  properties  of  orifice  flows  since  micronozzle  expansions  may  sufferAigh  viscous  losses  at  low 
pressure  operation.  The  mass  flow  and  relative  thrust  through  an{UDderexpanded)orifice  is  measured  as  a 
function  of  orifice  stagnation  pressure  from  0.1  to  3.5  Torr.  Nitrogen,  argon,  and  helium  propellant  gases 
are  passed  through  a  1.0  mm  diameter  orifice  with  a  wall  thickness  of  0.015  mm  .  Near-free  molecule, 
transitional  and  continuum  flow  regimes  are  studfied.  •  The  relative  thrust  is  determined  by  measuring  the 
displacement  of  a  novel  thrust  stand  designed  primarily  for  low  operating  pressure  propulsion  systems.  It  is 
shown  that  the  thrust  stand  deflection  is  a  function  of  the  facility  background  pressure,  and  corrections  are 
made  to  determine  the  deflection  for  a  zero  background  pressure  for  a  nitrogen  propellant 


Nomenclature 

a  -  speed  of  sound  (m/s  ec) 

A  -  area  (m2) 

c(y)  -  constant  dependent  on  ratio  of  specific 
heats 

c’  -  mean  molecular  thermal  speed  (m/sec) 
CD  -  discharge  coefficient 
d  -  diameter  (m) 

g  -  gravitational  constant  (=  9.8  m/sec2) 

Isp  -  specific  impulse  (sec) 
k  -  Boltzmann’s  constant  (=  1.38  x  10'23  J/K) 
Kn  -  Rnudsen  Number 
m  -  mass  (kg) 

M  -  mass  flow  (kg/sec) 
n  -  number  density  (m'3) 
p  -  pressure  (Pa) 

Re  -  Reynolds  Number 
rp  -  radius  of  penetration  (m) 
t  -  thickness  (m) 

T  -  temperature  (K) 


a  -  orifice  transmission  probability 
y  -  ratio  of  specific  heats 
X  -  mean  free  path  (m) 

H  -  viscosity  (Ns/m2) 
p  -  mass  density  (kg/m3) 

A  -  thrust  stand  deflection  (arb.  units) 

subscripts 

b  facility  background 
frn  -  free  molecule 
L  -  limit  (theoretical  maximum) 
meas  -  experimentally  measured 
o  -  stagnation  region 
p  -  plenum 

t  -  orifice  or  nozzle  throat  property 
theor  -  theoretical 


*  -  orifice  plane  (sonic)  region 


1.  Introduction 

In  recent  years,  micropropulsion  systems  have 
been  developed  to  address  the  need  for  highly 
mobile  micro-  and  nanospacecraft.  A  wide  array 
of  concepts  will  require  the  expansion  of 
propellant  gases  through  microscale  geometries 
(e.g.  micronozzles).  Because  of  the  volume  and 
power  restrictions  associated  with  storing  or 
producing  high  pressures  on  micospacecraft, 
many  micropropnlsion  systems  will  operate  at 
relatively  low  pressures  in  the  transitional  flow 
regime.1  The  Reynolds  number  gives  a  measure 
of  the  flow  efficiency  in  terms  of  viscous  losses. 
The  Reynolds  number  at  a  nozzle  throat  or  an 
orifice  is  given  by 

p,*  ,  LlA  m 


/tower  Reynolds  number  implies  higher  viscous 
flow  losses.  Microspacecraft  propulsion  systems 
may  inherently  operate  in  low  Reynolds  number 
regions  due  to  relatively  low  operating  pressures 
and  small  characteristic  dimensions. 


degrade  performance,  aid  expansion  of 
propellant  from  a  simple  thin  walled  orifice  may 
be  a  good  compromise  between  efficiency  and 
system  complexity.  Some  microelectro¬ 
mechanical  systems  (MEMS)^fabricated  nozzle 
geometries  involve  planar  or  rectangular  throats 
(i.e.  not  conical).3  Numerical  studies  have 
shown  that  flows  generated  near  the  side  walls 
can  result  in  even  higher  inefficiencies.6 

The  flow  complexities  from  sonic  orifices  have 
been  studied  for  several  years.7'11  However,  die 
determination  of  the  thrust  generated  from  gas 
expanding  through  an  orifice  is  an  area  that  has 
received  little  attention.  The  advent  of 
micropropulsion  systems  has  renewed  interest  in 
the  determination  of  propulsive  properties  of 
orifice  flows  since  micronozzle  expansions 
appear  to  have  major  viscous  losses. 


This  manuscript  explores  the  thrust  generated  by 
an  orifice  expansion  at  relatively  low  Reynolds 
numbed  Besides^  propulsion  system  in  its  own 
right,  these  orifices  are  also  being  investigated  as 
a  reliable  means  of  calibrating  micro-Newton 
thrust  stands.12 

2.  Theory 


S^tsJr 
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Figure  1  shows  the  specific  impulse  as  a  function 
of  distance  through  a  conical  micronozzle 
geometry  with  a  throat  diameter  of  d,  =  27.7  pm. 
Navier-Stokes  and  Direct  Simulation  Monte 
Carlo  numerical  simulations  have  been 
performed  at  two  different  stagnation  pressures 
i  ,  "*p0  =  106  Pa  and  10s  Par'*  For  cold  gas  operation 

i4  v  (T0  =  300  K).  the  corresponding  Reynolds 
numbers  are  1300  and  1*0  respectively.  As  Fig. 


1  shows,  the  specific  impulse  at  the  nozzle  exit  is 
approximately  14%  higher  than  at  die  nozzle 
throat  for  Re*  =  1300.  However,  the  specific 
impulse  at  the  exit  is  only  about  5%  higher  than 
at  the  nozzle  throat  for  Re*  =  130.  The  reduction 
of  efficiency  in  the  micronozzle  geometry  as  the 
Reynolds  number  decreases  is  due  to  viscous 
losses  near  the  nozzle  walls. 

There  is  currently  a  large  effort  being  devoted  to 
the  fabrication  of  micronozzles  with  throat 
diameters  on  the  order  of  one  to  tens  of 
micrometers  (microns).3'3  As  Fig.  1  indicates, 
nozzle  expansions  may  not  be  justified  when 
weighing  the  increase  in  performance  versus 
fabrication  complexity.  At  low  Reynolds 
numbers,  the  micronozzlc  geometry  may  in  fact 


To  assess  the  performance  of  the  orifice 
expansion  in  terms  of  propulsive  parameters, 
properties  at  the  entrance  plane  of  the  orifice  are 
calculated  from  known  stagnation  values.  The 
ratios  of  pressure,  density,  temperature  and 
velocity  for  inviscid  flow  are 
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The  theoretical,  inviscid  flow  value  for  the 
orifice  mass  flow  is 

M  =  p*a  At  ( 6) 
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Viscous  effects  can  be  measured  in  terms  of  a 
discharge  coefficient  defined  by 


(14) 
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where  Maili0r  is  calculated  from  Eq.  (6).  The 
theoretical  thrust  produced  bv  the  orifice  is  then 
given  by 

3  =  iln  +  p  At  (8) 
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(9) 


The  constant  c(v)  is  equal  to  1.16  and  1.14  for  y 
=  1.4  and  1.67,  respectively.  Therefore,  it  is 
expected  that  the  thrust  produced  by  the  orifice 
gas  flow  is  relatively  independent  of  the 
propellant.  ■■■  ■  ■ 


A  measure,  of  the  propulsive  efficiency  is  given 
by  the  specific  impulse  as 


!2(Y  +  l)k 
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For  free  molecule  fl-' .  the  Knudsen  number 
defined  by 

f-'!:  -  -f-  (11) 

dt 


is  relatively  high  (Kn  >  10).  This  is 
accomplished  at  very  low  stagnation  pressure 
operation  where  the  molecule  mean  free  path  is 
larger  than  the  .orifv;  diameter.  The  free 
molecule  mass  flow,  thrust  and  specific  impulse 
are  given  by 


(12) 


(13) 


3.  Experiment 

The  orifice  used  in  this  study  has  a  diameter  of  d, 
=  1.0  mm  and  a  thickness  of  t,  =  0.015  mm 
giving  a  t/d  =  0.015.  For  t/d  =  0.015,  the 
transmission  probability  (a  in  Eqs.  (12)  and  (13)) 
is  very  close  to  unity.13  The  orifice  is  machined 
by  conventional  means  in  a  Tantalum  shim 
which  is  attached  to  an  aluminum  plenum  as 
shown  in  Fig.  2.  The  aluminum  plenums  are 
attached  to  a  torsional  thrust  stand  shown  in  Fig 
3.  The  thrust  measurements  involve  sensing  the 
angular  displacement  resulting  from  a  torque 
(thrust  force)  applied  to  a  damped  rotary  system. 
The  present  method  for  detecting  angular 
deflection  is  to  measure  the  linear  displacement 
of  a  known  radial  distance  nsing  a  linear  voltage 
differential  transducer  (LVDT)  by  Macro 
Sensors.  The  total  linear  movement  of  the  arm  is 
approximately  0.5  mm  for  a  2mN  thrust  level 
which  corresponds  to  less  than  0.1°  angular 
deflection.  The  detailed  operational 
characteristics  of  this  thrust  stand  is  the  topic  Of 
earlier  work.7 

The  thrust  stand  is  placed  inside  the  CHAFF-II 
facility,  a  steel  vacuum  chamber  pumped  by  a 
Roots  blower  system  with  a  pumping  speed  of 
2000  L/sec  for  nitrogen.  Ultimate  pressures 
achievable  in  CHAFF-II  are  approximately  1.0  x 
10-4  Torn. 

The  propellant  is  introduced  into  the  orifice 
plenum  through  an  adjustable  needle  valve 
located  downstream  of  an  MKS®  mass  flow 
meter.  In  the  experimental  configuration,  the 
mass  flow  meter  operated  in  the  continuum 
regime  through  the  pressure  range  studied. 
Nitrogen,  argon  and  helium  are  used  as 
propellant  gases  in  this  study. 

4.  Results 

.  j  • 

Figure  4  shows  the  discharge  coefficient  as  a 

function  of  the  Reynolds  number  for  nitrogen. 

At  lower  Reynolds  numbeff  the  measured  values  $ 
should  asymptotically  approach  the  theoretical 
free  molecule  limit  of  0.583.  For  die  range  of 
Reynolds  numbed  shown  in  Fig.  4,  die  orifice  $ 
stagnation  pressure  ranges  from  0.1  to  3.5  Tore. 

As  expected,  the  discharge  coefficient 
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asymptotes  towards  unity  for  higher  Reynolds 
number 
a 

The  measured  linear  deflection  from  the  thrust 
stand  mounted  LVD'f  is  shown  as  a  function  of 
orifice  stagnation  pressure  for  nitrogen 
propellant  in  Fig.  5.  Using  Eq.  (8)  with  the 
measured  mass  flow  and  stagnation  pressure,  the 
theoretical  range  of  thrust  shown  in  Fig.  5  is 
from  approximately  7.9  /iN  (p0  =  0.1  Torr)  to 
430  juN  (p0  =  3.5  Torr).  Similar  plots  are  shown 
for  argon  and  helium  in  Figs.  6  and  7, 
respectively.  Figure  8  shows  the  nitrogen 
deflection  as  a  function  of  mass  flow. 

As  seen  in  Fig.  9a,  the  data  for  nitrogen  and 
argon  can  be  fit  by  the  same  straight  line  as 
expected  from  Eq.  (9).  However,  this  level  of 
agreement  is  only  true  for  propellants  which 
have  a  similar  Reynolds  number  for  a  given 
stagnation  pressure  (i.e.  similar  viscous  effects  as 
a  function  of  p„).  Because  the  orifice  Reynolds 
number  (Eq.  (1))  is  a  factor  of  three  lower  for 
helium  than  nitrogen  at  a  given  stagnation 
pressure,  it  is  expecied  that  the  thrust  (deflection) 
would  be  somewhat  lower  for  helium'  due  to 
viscous  effects  as  shown  in  Fig.  9b. 

The  data  for  helium  and  nitrogen  are  shown  in 
Fig.  10.  The  discrepancy  in  the  data  is  indicative 
of  the  fact  that  the  helium  How  is  more  rarefied 
(higher  Knudsen  number)  for  a  given  stagnation 
pressure.  For  example,  for  nitrogen  flow  at  p0  = 
0.1  Torsuhe  Knudsen  number  is  approximately 
0.5.  For  helium,  the  Knudsen  number  is 
approximately  1.4  for  the  same  stagnation 
pressure.  This  indicates  that  viscous  effects  are 
more  important  for  helium  flow  than  nitrogen  at 
a  given  stagnation  pressure. 

It  is  known  that  the  deflection  for  a  given  orifice 
stagnation  pressure  is  dependent  on  the 
background  pressure  <>:'  ;hc  facility.  Figure  11 
shows  the  measured  de!  lection  for  given  orifice 
stagnation  pressures  as  a  function  of  the  chamber 
background  pressure.  The  absolute  deflection 
for  a  given  stagnation  pressure  decreases  as  the 
background  pressure  increases.  The  slope  and 
intercept  of  the  data  is  used  to  correct  the  data  in 
Fig.  5  to  a  zero  background  pressure  in  the 
following  section.  • 

5.  Discussion 


Indications  of  Rarefied  Flow 

Figure  12  shows  the  effects  of  rarefied  flow  on 
the  measured  thrust  (deflection)  from  the  orifice. 
For  helium,  the  flow  Knudsen  number  at  p0  =  0.1 
Ton-  is  approximately  1.4.  The  helium  deflection 
data  in  Fig.  12  shows  the  theoretical  lines  for 
free  molecule  and  inviscid  continuum  thrust 
The  data  is  bounded  by  the  free  molecule  and 
continuum  solutions.  At  the  lower  operating 
pressure,  the  data  closely  follows  the  free 
molecular  slope  and  tends  toward  the  inviscid 
continuum  solution  at  higher  pressures.  The  data 
does  not  quite  asymptote  to  the  inviscid  solution 
at  the  maximum  Reynolds  number  (p0  =  3.5 
Torr)  of  about  27. 

Figure  13  shows  a  similar  result  for  nitrogen 
propellant  flows.  At  the  higher  operating 
pressures,  the  nitrogen  flows  closely  follow  the 
inviscid  continuum  solution:  This  is  consistent 
with  tile  discharge  coefficient  results  shown  in 
Fig.  4  which  indicates  a  discharge  coefficient 
near  unity  for  Re*  >  40. 

Facility  Background  Pressure  Corrections 

Because  the  facility  background  pressure  is  made 
up  of  two  components  (laboratory  air  and 
propellant),  corrections  to  the  thrust  stand 
deflections  can  only  be  approximated  for  the 
nitrogen  propellant  cases  where  the  propellant 
and  the  laboratory  backgrounds  are  similar.  * 

The  mechanism  for  lower  thrust  deflection  as  a 
function  of  increased  background  pressure  is 
shown  schematically  in  Fig.  14.  For  no  orifice 
flow,  the  background  pressure  exerts  an  equal 
force  on  the  front  and  back  sides  of  the  orifice 
plenum  (equilibrium).  As  flow  is  introduced 
through  the  orifice,  the  resulting  jets  acts  like  an 
ejector  pump  similar  to  that  of  the  oil  vapor  in  a 
vacuum  diffusion  pump.13-14  Collisional 
“removal”  of  the  background  gas  by  the  orifice 
plume  results  in  a  lower  background  pressure  on  ^ 
the  jet  side  of  the  plenunme5  on  the  back  side. 
This  pressure  difference  exerts  a  force  on  the 
orifice  in  a  direction  opposite  of  the  thrust  vector 
prodneed  by  the  jet.  Since  the  gas  density  in  the 
plume  is  relatively  high  compared  to  the 
background  gas  density  in  the  vicinity  of  the 
orifice,  the  source  flow  can  effectively  prevent 
background  molecules  from  penetrating  the 
orifice  plume  and  striking  the  front  surface  of  the 
orifice  plenum  for  reasonable  background 
pressures  (pb  <  1  x  10’3  Torn).15  This  suggests 
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that  the  deflection  as  a  function  of  the 
background  pressure  should  be  linear  as  a  first 
order  approximation  as  Fig.  1 1  indicates,  or 

dA 

- =  constant 

fyb 

{^Ptso)~[^pt)  =  Pb^ejf  ^ 

The  effective  area,  A^,-,  is  a  measure  of  the  area 
being  utilized  to  cause  the  deflection  opposite  to 
the  thrust  vector.  Assuming  that  the  background 
pressure  in  front  of  the  orifice  plenum  is  zero, 
the  force  exerted  on  the  back  side  of  the  orifice 
plenum  (Ap  =  22.75  cm’)  is  approximately  100 
/tN  for  p„  =  3.3  x  10''  Tour.  Therefore,  the 
background  pressure  “negative  thrust”  effect  can 
be  a  significant  fraction  of  total  measured  thrust 
produced  by  the  orifice. 

Figure  15  shows  the  nitrogen  thrust  stand 
deflection  for  a  corrected  zero  background 
pressure.  This  data  is  derived  from  the  slopes  of 
the  deflection  versus  background  pressure  curves 
in  Fig.  11  and  the  data  obtained  in  Fig.  5. 

6.  Conclusions 


7.  Acknowledgements 

The  authors  wish  to  thank  Mr.  Brian  Bjelde  for 
his  assistance  with  the  thrust  stand  and  data 
reduction. 

8.  References 

1.  Ketsdever,  A.,  “System  Considerations  and 
Design  Options  for  Microspacecraft  Propulsion 
Systems,”  in  Micronroptilsion  for  Small 
Spacecraft.  AIAA  Progress  Series  in 
Astronautics  and  Aeronautics,  eds.  Micci  and 
Ketsdever,  Vo!.  187, 2000,  pp.  139-166. 

2.  Ivanov,  M.,  Markelov,  G.,  Ketsdever,  A., 
Wadsworth,  D.,  “Numerical  Study  of  Cold  Gas 
Micronozzie  Flows,”  AIAA  paper  99-0166, 
January,  1999. 

3.  Bayt,  R.,  Breuer,  K.,  “Fabrication  and  Testing 
of  Micron-Sized  CuKi-Gus  Thrusters,”  in 
Micropromilsion  for  Small  Spacecraft.  AIAA 
Progress  Series  in  Astronautics  and  Aeronautics, 


eds.  Micci  and  Ketsdever,  Vol.  187,  2000,  pp. 
381-398. 

4.  Janson,  S.,  Helvajian,  H.,  “Batch-Fabricated 
Microthrusters:  Initial  Results,”  AIAA  paper  96- 
2988,  July  1996. 

5.  Kohler,  J.,  Jonsson,  M.,Simu,  U.,  Stemnark, 
L.,  “Mass  Flows  of  Fluidic  Components  for  Cold 
Gas  Thruster  Systems,”  Micro/Nanotechnology 
for  Space  Applications,  Pasadena,  CA,  April 

1999. 

6.  Alexeenko,  A.,  Gimelshein,  S-,  Levin,  D., 
Collins,  R.,  “Numerical  Modeling  of 
Axisymmetric  and  Three-Dimensional  Hows  in 
MEMS  Nozzles,”  AIAA  paper  00-3668,  July 

2000. 

7.  Ashkenas,  H.,  Sherman,  F.,  “The  Structure 
and  Utilization  of  Supersonic  Free  Jets  in  Low 
Density  Wind  Tunnels,”  in  Rarefied  Gas 
Dynamics.  Proceedings  of  the  4th  International 
Symposium,  ed.  de  Leeuw,  Vol.  2,  1966,  pp.  84- 
105. 

8.  Sreekanth,  A.,  Prasad,  A.,  Prasad,  D., 
“Numerical  and  Experimental  Investigations  of 
Rarefied  Gas  Flows  Through  Nozzles  and 
Composite  Systems,”  in  Rarefied  Gas  Dynamic?, 
Proceedings  of  the  IT®  International  Symposium, 
ed.  Beylich,  1991,  pp.  987-994. 

9.  Rebrov,  A.,  “Free  Jets  in  Vacuum 
Technologies,”  J.  Vac.  Sci.  and  Technol.  A,  to  be 
published. 

10.  Meyer,  J.,  “Measurements  of  Particle 
Densities  and  How  Directions  in  Freejets  and  in 
Background  Gas,”  in  Rarefied  Gas  Dynamics: 
Experimental  Techniques  and  Physical  Systems. 
AIAA  Progress  Series  in  Astronautics  and 
Aeronautics,  eds.  Shizgal  and  Weaver,  Vol.  158, 
1994,  pp.  333-341. 

11.  Livesey,  R.,  “How  of  Gases  Through  Tubes 
and  Orifices,”  in  Foundations  of  Vacuum 
Science  and  Technology,  ed.  Lafferty,  1998,  pp. 
81-140. 

12.  Tew,  J.,  VanDenDriessche,  J.,  Lutfy,  F., 
Muntz,  EJ>.,  Wong,  J.,  Ketsdever,  A.,  “A  Thrust 
Stand  Designed  for  Performance  Measurements 
of  the  Free  Molecule  Micro-Resistojet,”  AIAA 
paper  2000-3673,  July  2000. 


5 


13.  Dayton,  B.,  “Diffusion  and  Diffusion- 
Ejector  Pumps,”  in  Foundations  of  Vacuum 
Science  and  Technology,  ed.  Laffertv,  1998,  pp. 
175-232. 

14.  Rebrov,  A.,  “Studies  on  Physical  Gas 
Dynamics  of  Jets  as  Applied  to  Vacuum  Pumps,” 
in  Rarefied  Gas  Dynamics.  Proceedings  of  the 
15t''  International  Symposium,  ed.  Teubner, 
1986,  pp.  455-473. 

15.  Muntz,  E.P.,  Hamel,  B.B.,  Maguire,  B.L., 
“Some  Characteristics  of  Exhaust  Plume 
Rarefaction,”  AIAA  J.,  Vol.  8,  No.  9,  1970,  pp. 
1651-1658. 


viewport 


Pressure  (Torr) 


— s —  Deflection  11-15 
— a  -Deflection  11-16 
-  e-  -  Deflection  11-17 


>—■ . —  ———I 

1  _ 

2  2 
sure  (Torr) 


FM  Thrust  (N) 


- Normalized  FM  Thrust. 

-  -  Normalized  Continuum  Thrust 
o  Deflection  11-21 
x  Deflection 
+  •  -  Deflection  1 1-27  •  * 


'  X -p  pe&iiUjc, 
CU*c^  ULp-koO  . 


FM/Cont  d 


- y  =  8j0476e-11  +  D.045541x  R=  1 

- y  =  -i.7208e-10-4lo.11549x  R=  1  / 


Stagnation  Pressure  (Torr) 


Deflection 


L&jen, 


Jt  ? 


0.35 


0.30 


0.25 


0.20 


0.15 


0.10  h 


0.050  h 


- y  =  0.050114  -  37.089X  R=  0.9984 

I - -y  =  0.1 3448  -  91 ,42x  R=  0.99882 


— t  j  i  i  i  i  |  i  i  i  i  |  i  r 

y  =  0.01 4136-1 0.1 36x  Rr  0.98927  ^ 


-x 


X. 


-  y  =  0.2866  -  194.43x  R=  0.99863 
y  =  0.45203  -  278.07x  ^  0.99904 


A 


O-  _ 


-OO 


□ 


E3- 


■  -O- . . 


-1 _ I _ I _ 1—1 _ I _ ! _ I _ L 


_! _ _ L_ 


J _ I _ I _ I _ I _ I _ I _ 1 - ! - 1 - 1_ 


0.0 

0.0001  0.0002  0.0003  0.0004  0.0005  0.0006 


Row  Pressure  (torr) 


el  - 1 

Ff)-« 


0.00 


Normalized  FM  Thrust 


- Normalized  FM  Thrust 

—  -  Normalized  Continuum  Thrust 
o  Deflection  11-21 
x  Deflection  .11-20 


: - y  =  8i0476e-11 +p.045541x  R=  1  /j 

- V  =  -1.7208e-1 0  4  0. 1 1 549x  R=  1  / . 


\x 
. \ 

■/*’ 

* 

2  3 

Stagnation  Pressure  (T orr) 


as 


